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An Experimental Study Into The Dynamic 
Behavior of Hemispherical Shells Subjected 
To External And Internal Loading 

by 

Panagiotis Romanos 
Abstract 

An experimental investigation into the dynamic 
behavior of fully clamped hemispherical shells subjected 
to uniformly distributed external and internal impulsive 
loading is reported. The applied loads were sufficient 
to cause plastic flow of the material. The materials 
used were mild steel 1020 and aluminum 606I-T6 . The 
deflections of internally loaded specimens showed a lin- 
ear variation with the initial velocity, while the de- 
flections of the externally loaded specimens vary non- 
linearly with velocity after a certain value of the 
latter; this may be useful in evaluating a "buckling 
threshold" . The externally loaded specimens exhibited 
a number of modes of deformation in the meridional 
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direction which is much less than the number of modes 
of deformation in the circumferential direction. 
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INTRODUCTION 



The considerable attention which has been given 
in recent years to the problem of instability of shell 
structures, is primarily due to their increasing use 
in hydrospace vehicles, nuclear reactors and space ve- 
hicles. For example, complete spheres are used for 
deep submersibles, while spherical caps or hemispheres 
are used in order to form the ends of the pressure hulls 
of submarines. 

These hydrospace vehicles as well as other struc- 
tures are often subjected to high intensity loading which 
may be applied either dynamically or statically. Such 
high loadings can be caused by deep diving, collision, 
water wave impacts or high explosive effects, the latter 
being usually associated with naval vessels in wartime. 

Although considerable work has been so far devoted 
to instability of shells subjected to static loadings, 
the work done on dynamic loading is limited. For the 
case of hemispherical shells which is the subject of 
this investigation, a literature survey showed that the 
subject of instability of hemispherical shells subjected 
to external dynamic loading has been barely dealt with 
up to now. In Reference (1) some experimental work has 
been done on the dynamic behavior of hemispherical shells 
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subjected to dynamic external loading; two types of 
dynamic tests were examined, of which one involved blast - 
induced dynamic buckling and permanent deformation of the 
shell, while the second involved producing permanent de- 
formation by projectile impact coincident with the axis 
of symmetry of the hemispherical shell. Hence in both 
cases, the loading was not uniformly distributed over 
the surface of the shell. 

In Reference (2), theoretical work has been done 
on elastic-plastic buckling of shells under dynamic load- 
ing, with the case of hemispherical shells included. 

This is the only theoretical work on dynamic instability 
of hemispherical shells subjected to external loading 
that the author is aware of. 

It is the purpose of a major part of the present 
investigation to study the experimental behavior of 
clamped hemispherical shells subjected to externally 
applied impulsive loading. It is hoped that the results 
presented herein will help in developing new theories 
on this subject or confirming already existing theories 
as those of Reference (2). 

The remainder of this work is devoted in the in- 
vestigation of the dynamic behavior of clamped hemisphe- 
rical shells subjected to internal loading. This forms 
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f 



a continuation of work reported in References (5), (4) 
and (5) and differs in the specimen geometry (primarily 
the ratio of radius to the thickness of the hemisphere) . 
The dynamic internal loading of hemispherical shells has 
an important application in the case of an internal ex- 
plosion occurring in pressure vessels or reactor vessels. 



- 3 - 



Notation Used 



nominal explosive thickness 
h measured average explosive thickness 

€L 

H average thickness of hemispherical specimen 
I g p specific impulse of the explosive 

p 

Mqs equal to a QS H /4 , is the static collapse 
moment of the hemispherical specimen 

m g loaded mass of the hemispherical specimen 

R g mean radius of the hemispherical specimen 

R g radius formed between the hemispherical part 
of the specimen and the flange 

V QS impulse velocity imparted to hemisphere 

W weight of explosive used 

V 

W permanent radial deflection, maximum value of 
deflections for internally loaded specimens 
and average value of deflections for externally 
loaded specimens 

p density of specimen material 

* s 

a average static yield stress of hemispherical 

o s 

specimen 
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Notation Used (cont’d) 



a (V ) 2 R 2 

equal to ■■ os ^ ' OS , is a non-dimensional 



os 



parameter used to relate results to one another 



2a overall angle of the hemisphere which is covered 
by explosive 
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EXPERIMENTAL PROCEDURE 



Twenty eight hemispherical specimens were used in 
this study. The general dimensions of these specimens are 
shown in figure 1 . Twenty one of the specimens were made 
from hot rolled mild steel 1020 and the remaining seven 
from aluminium 6061-T6. The specific dimensions for each 
specimen used are shown in table 1. The specimens were 
initially hydroformed from plates of 1/4" thickness, using 
a Cincinatti Milling Hydroform machine, which gave the 
hemispherical shape to the flat plates. The aluminium 
specimens were hydroformed while being in the zero heat 
treatment condition and were heat treated after the 
hydroforming, the T6 condition . The steel specimens were 
annealed after hydroforming. Then all the specimens were 
machined in order to be brought down to the dimensions 
required. Unfortunately, although the initial order was to 
machine the specimens at a maximum thickness of 0.145", due 
to an error of the machining company the maximum thickness 
obtained was only 0.090". Thus the specimens delivered 
were 10 steel specimens of 0.070" thickness, 11 steel 
specimens of 0.090" and 7 aluminium specimens of 0.090" 
thickness . 

Initially the specimens were measured by the 
apparatus shown in figure 7. This apparatus consists of a 
dial gage with stem travel of 3" and with subdivisions of 
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0.001", which is attached permanently on the head of a 

milling machine. Each specimen was marked before the 

measurement, so that in the measurements taken after 

deformation points corresponding to the same meridional 

direction could be measured. The axis of the dial gage 

was inclined by 20° to the vertical direction, so that 

the part of the hemisphere which is near the flage can be 

measured. The table on which the specimen to be measured 

was resting could be rotated 360° in increments of one 

degree and it could travel longitudinally with increments 

of 0.001". The specimen was placed in the center of the 

table and the dial gage was arranged in such a way that its 

needle would pass through the center of the revolving plate 
and therefore through the top of the hemispherical specimen 

when it was centered on the plate. The specimens were held 

firmly on the table by three cap screws placed every 90 

degrees. For the specimens which were to be externally 

loaded initial measurements were taken every 30° degrees, 

and in each such position of the specimen 20 points were 

measured starting from a point nearest to the edge of the 

flange and ending to the top of the hemispherical part. 

For the specimens which were to be internally loaded 

initial measurements were taken every 45 degrees and in 15 

points in the meridional direction up to the top of the 
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hemispherical part. The dial gage was adjusted to indicate 
zero on the top of the plate where the specimen was 
resting, so that the measurements gave the distance from 
the surface of the specimen to the top of the plate. In 
order to get the vertical distance of each point measured 
from the top of the plate the dial gage indication was 
multiplied by cos 20° since as mentioned above the dial 
gage had an inclination of 20° to the vertical . Similarly , 
the longitudinal position of each measured point was found 
by subtracting the product of the dial gage indication 
times sin 20° from the longitudinal travel of the table 
which was measured from the edge of the specimen's flange. 

A further explanation of how these measurement conversions 
were made is given in appendix I. The thickness of the 
specimen was measured by the device shown in two views in 
figures 8 and 9. Thirty eight measurements were taken in 
each specimen, in two intersecting at right angles meridional 
directions. Of the nineteen measurements in each of these 
directions eleven were on the hemispherical part and the 
remainder 8 in the flange. The twenty two measurements 
taken on the hemispherical part of each specimen were 
averaged and this average thickness is shown as H in table 
1. From this table it can be seen that only specimens 13 
17 and 18 have an average thickness which is slightly less 
than the required, which would be 0.088 since the required 



- 8 - 



thickness tolerance was 0.002" and the nominal thickness 
of the mentioned specimens was 0.090". 

After the initial measurements were taken, the 
area of the specimen which were to be exposed to the 
dynamic load was covered by polyurethene foam, of 1/4 inch 
thickness, of polyester type, with density 0.032 grams/cm3. 
This foam was glued in place with Dupont # 4684 cement, and 
it acted as an attenuator which prevents spalling of the 
material which would happen due to the explosion. Next, the 
explosive was cut in four pieces which had the form of 
spherical triangles and the thickness of each of the pieces 
cut was measured on seven points (.six around the perimeter 
and one in the middle) for tests 1 to 24 and on four points 
(three in the corners and one in the middle) for tests 25 
to 28. The device with which the explosive thickness was 
measured is shown in figure 10 and it consists of a flat 
plate made of glass and a dial gage placed with its stem 
vertical to the plate and with its measuring point 
substituted by a wooden pad, so that it will not exert a 
localized high pressure on the explosive. The explosive 
thickness were found varying from 0.0060" to 0.0125" for 
explosive of nominal thickness 0.010", from 0.0152" to 
0.0196" for explosive of nominal thickness 0.015" and from 
0.028" to 0.0335" for explosive of nominal thickness 0.030". 
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The average thickness of the explosive used in each test 
is also shown in table 1. The explosive was then weighed 
with a scale of accuracy 0.01 gram and was glued on the top 
of the foam attenuator. In all tests from 1 to 24 the 
edge of the explosive after it was glued on the foam had 
a vertical distance from the top of the flange of 1/4" for 
the externally loaded specimens and 1/4" from the bottom 
of the flange for the internally loaded specimens; from this 
distance the angle formed by the radii of the spherical 
surface which pass from the edges of the explosive was 
found for each specimen after drawing the initial profile 
of the specimen as described in appendix I. The angles so 
measured are shown in table 1. In tests 25 to 28 the angle 
mentioned above was specified and is also shown in table 1. 

After the explosive had been firmly glued on the 
foam attenuator, the specimen was bolted in the clamping 
devices which are shown in figures 2 and 3. Figure 2 shows 
the clamping device used for external loading and figure 3 
shows the clamping device used for internal loading. Both 
clamping devices were case hardened in order to protect 
the serrations which prevented the slippage of the specimens. 
The specimen was bolted between the two parts of the 
clamping devices with 16 heat treated steel bolts of 1/4" 
diameter and twenty threads per inch. The lower part of the 
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clamping device was then secured on the test table which 
is shown in figure 4 by 4 small steel clamps as shown in the 
arrangement of figure 5. Figure 5 shows the arrangement 
used for internal loading. The only difference in the 
external loading arrangement is that the hemispherical part 
of the specimen is upwards instead of downwards as shown 
in figure 5. 

The explosive used in all the experiments came 
from the same order and it was calibrated as described in 
detail in appendix II. An explosive leader of dimensions 
30 inches by 0.25 inches by 0.015 inches was used in the 
experiments. One end of the leader over an area of 0.5"X0.25" 
was pressed firmly on the middle of the hemispherical 
specimen in such a way as to join the four parts of the 
explosive which had been previously glued on the foam 
attenuator. The other end of the explosive leader was folded 
over a length of 1.5 inches and on this end was pressed a 
number six blasting cap which detonated the the explosive. 
The blasting cap was put into the hole of a metal block under 
which the end of the explosive leader was placed. With this 
arrangement the specimen was not affected by the impulse due 
to the blasting cap. The overall testing arrangement is shown 
in figure 6 with the specimen arranged for external loading. 

After each experiment, the deformed specimen was 
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measured with the same apparatus which was used for the 
initial measurements as mentioned above and is shown in 
figure 7. Since the external surface of the specimens were 
used for the measurements, the externally loaded specimens 
were washed with pure aceton, so that all remaining parts 
of burned attenuator were removed and the specimens were 
completely clean during the measurements. In these 
measurements only the points along one semicircle were 
measured, for which the initial measurements had shown the 
greatest variation between the two parts of the semicircle; 
this corresponds to the worst possible variation in the 
radius measurements made. For the internally loaded specimens 
the same number of points were measured after deformation 
as the ones of the initial measurements. For the externally 
loaded specimens, because of the buckling phenomenon, 
irregular shapes were observed and hence eighty points were 
measured on a semicircle, so that an accurate deformation 
profile could be drawn. In addition, for the externally 
loaded specimens measurements were taken at two distances 
from the edge of the flange, every 10 degrees around the 
periphery so that deformation profiles in the sense of circles 
which before deformation were parallel to the base could 
be obtained. The details of how the deformation profiles were 
obtained and drawings of deformation profiles for each of the 
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specimens are given in appendix I . 

The data obtained with this experimental 
investigation are shown in tables 1, 2 and 3 and plotted 
in figures 11, 12, 13, and 14. Sample calculations for 
obtaining these data are shown in appendix I. 
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FIGURE 1 : General specimen geometry 
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FIGURE 2: Clamping device for external loading 
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FIGURE 3. Clamping device for internal loading 
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Figure 8 
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Figure 9 
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Figure 10 
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Remarks On TABLE 2 




Test 

No. 

2 

4 

6 

11 

12 

13 

14 

19 

21 



TABLE 3 



MODES OF DEFORMATION OF EXTERNALLY 
LOADED SPECIMENS 



No. of Peaks 
In Meridional 
Semicircle 



No. of Peaks 
In Circles Parallel 
To The Base 



3 



11 



5 

3 

3 

3 

4 



3 

3 



10 

7 

10 

8 

11 

8 

9 

10 
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Experimental Results 



The results obtained in the present experimental 
investigation are shown in Tables 2 and while some 
experimental data used to evaluate those results are 
given in Table 1. For the calculation of the quantities 
shown in Table 2 the yield stress of the material used 
for the fabrication of the specimens, as well as their 
densities were required; for this purpose tensile tests 
were conducted with tensile specimens from the materials 
used and are reported in Appendix III. For the densities 
of the materials average values of densities of the same 
materials obtained in previous investigations were used, 
and are also reported in Appendix III. Samples of the 
calculations made in order to obtain the results shown 
in Table 2 are shown in Appendix I. In Appendix I the 
method of measuring the radial deflections is described, 
tables of these deflections for each specimen are given 
and deformation profiles in the meridional direction for 
one measured semicircle are given for each specimen. 

From the deformation profiles of the externally loaded 
specimens in the meridional direction, the number of peaks 
of the deformation profile were measured and are shown in 
Table 3 . For the externally loaded specimens two deform- 
ation profiles in the circumferential sense were also 
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obtained for each specimen by the method described in 
Appendix I and they are drawn also in the same Appendix 
for each of the externally loaded specimens. The last 
mentioned circumferential profiles were also drawn for 
specimen 14 although the radial deformations of the speci- 
men were not measured because of a large crack which 
developed in the part between the spherical part and the 
flange of the specimen; these profiles were drawn be- 
cause it was considered useful to see the deformation pro- 
files in the circumferential direction of this specimen 
which was subjected to higher loading than the other speci- 
mens of the same nominal thickness. 
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Discussion of Results 



The deflections of the specimens used in the meri- 
dional and circumferential direction which are presented 
in form of tables and drawings in Appendix I show a cer- 
tain non symmetry. This may be due to a number of factors. 
The radius of the hemisphere measured on one semicircle in 
the meridional direction was found varying a little so that 
the shape was not perfectly circular. The thickness was 
found varying by a little more than +0.002" which was the 
tolerance initially required; moreover the thickness of the 
specimens in the region between the hemispherical part and 
the flange was less than the thickness of the rest of the 
specimen and this can justify the comparatively high de- 
flections measured in this region, since one would expect 
to have deflections increasing as we go away from the 
clamped region and not decreasing near this region and then 
increasing again. The main factor to which the non sym- 
metry of deformation may be due to the variation in the 
thickness of the explosive which is mentioned in the Ex- 
perimental Procedure. 

Since the energy ratios E r calculated for all 
specimens in the way shown in Appendix I and reported 
in Table 2 vary from 43 to 425, the material elasticity 
should not affect the overall structural response and the 
energy dissipation can be considered as purely plastic. 

-4o- 



From the deformation profiles for the externally- 
loaded specimens which are presented in Appendix I, the 
conclusion may be drawn that the number of modes of de- 
formation in the meridional direction is much less than 
in the circumferential direction as shown also in Table 3. 

The experimental results are plotted in Figures 11 to 
l4. In Figure 11 the deflection of the specimen is plot- 
ted versus the velocity imparted to the specimen by the 
explosive. For the externally loaded specimens the aver- 
age deflection was used while for the internally loaded 
specimens the average deflection was used. It can be seen 
from Figure 11 that for the same velocity V , the speci- 

OS 

mens with nominal thickness 0.070 in. exhibited deflection 
higher than those of nominal thickness 0.090 in., which, 
of course, was to be expected. Moreover, although the de- 
flections of the internally loaded specimens show an almost 
linear variation with velocity, the deflection of the ex- 
ternally loaded specimens shows a strong non-linear varia- 
tion with velocity after a certain point. This non-linear- 
ity is more pronounced in Figure 12 in which the same plot- 
ting as in Figure 11 is made, but using the maximum deflec- 
tions of the externally loaded specimens. This implies that 
a "buckling threshold" could be found, which is the point 
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after which the slope of the deflection-velocity curve 

starts to change abruptly. 

In Figure 13 the non dimentional deflection W/H 

is plotted against the non-dimensional impulse parameter 

w 

X . The values of -jp and X for each specimen were 

taken from Table 2. It can be seen from Figure 13 that 

W 

the variation of the non-dimensional deflection -^p with 
the non-dimensional impulse parameter X is slightly 
non-linear. Also for the same value of X the values 
of -p- for the specimens with nominal thickness of 
0.090 in., are higher than those for the specimens with 
a nominal thickness of 0.070 in. 

Finally in Figure 14 the non-dimensional deflection 
W 

-jp is plotted against half of the included angle of the 
explosive a , for specimens 26 and 28, which were the 
only aluminum specimens for which results were obtained. 
These experimental points are too few for valid conclu- 
sions to be drawn. In Figure 14 the theoretical pre- 
dictions for the values of which were calculated in 

n 

Appendix I according to a method given in Reference (9) 
are also plotted. It can be seen that the theoretical 
prediction for the value of -p- for specimen 26 is much 
less than the experimental value, while for specimen 28 
a much better agreement was obtained. 
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CONCLUSIONS AND RECOMMENDATIONS 



An experimental investigation into the dynamic 

behavior of fully clamped hemispherical shells subjected to 
uniformly distributed external and internal impulsive loading 

is reported herein. The applied loads were sufficiently high 

to cause plastic flow of the material. The materials used 

for the specimen fabrication were mild steel 1020 which is 

a highly strain rate sensitive material and aluminium 

6061-T6 which is insensitive to strain rate. It was one of 

the purposes of this work to assess the influence of strain 

rate sensitivity of the materials on the structural response. 

However due to the small thicknesses of the delivered 

specimens this became impossible because complete failure 

of the aluminium specimens was observed at comparatively 

low loads. 

It was observed that the deflections of the 
internally loaded specimens show a linear variation with 
the initial velocity imparted to the specimens by the 
explosive, while the deflections of the externally loaded 
specimens vary non-linearly with velocity after a certain 
value of the latter. This may imply that a "buckling 
threshold" could be calculated as the velocity for which 
the slope of the deflection vs velocity curve starts to 
change abruptly. 

From the results of externally loaded specimens 
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it was observed that the number of modes of deformation in 
the meridional direction is much less than that in the 
circumferential direction. 

It is recommended that some more tests be made 

with steel specimens of 0.070 in thickness so that the 

results obtained in this work can be verified, because for 

this thickness the number of successful tests made were 

only a few. Also it would be interesting to test thicker 

aluminium specimens with external load as well as steel 
specimens with the same thickness in order to find the 

behavior of aluminium specimens in buckling and to assess 

the influence of strain rate sensitivity of the material 

on the deformation; specimens of higher thicknesses will 

permit also the use of higher loadings without complete 

failure . 
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APPENDIX I 

EXPLANATION OP MEASUREMENT CONVERSIONS, DRAWING OF 
INITIAL AND DEFORMATION PROFILES TABLES OF DEFLECTIONS 
AND SAMPLE CALCULATIONS. 
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Explanation of Measurement Conversions 



The manner and the apparatus with which the speci- 
mens were measured before and after deformation were 
described above. Figure 1-1 shows a schematic diagram 
of the measuring apparatus. In this figure the solid 
arrow on the spherical surface indicates the stem of 
the dial gauge with which the measurements were taken 
and which had an inclination of 2C ? to the vertical as 
shown. The dial gauge reading is the distance d , 
that is the distance from the spherical surface to the 
top of the flat plate on which the specimen was resting. 
The longitudinal travel of the table on which the flat 
plate and the specimen were resting is shown as l in 
Figure 1-1 and as mentioned before it could be measured 
with accuracy of 0.001". In order to plot the initial 
and deformation profiles, the distances x and y which 
are shown in Figure 1-1 were required. Distance x is 
the horizontal distance of the measured point from a point 
near the edge of the flange which was the zero of the 
longitudinal travel of the table. From the geometry of 
Figure 1-1 we get: 

y = d cos 20° 

x = l - d sin 20° 
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For example with d = I.158I in. and l = 1.875 in. 
which were taken from initial measurements of specimen 
11, we have: 

y = 1.0082 in. and x = 1.4789 in. 

In this way all the measurements taken before and 
after deformation were converted and the coordinates of 
each point obtained, so that they could be easily plotted. 
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Drawing of Initial and Deformation Profiles 



After the measurements taken before and after de- 
formation had been converted as described above, the 
initial and deformation profiles of one semicircle (of 
which the before deformation measurements showed the 
highest variation) were drawn with a scale of 10 to 1 . 
From the initial profile so drawn the outside radius of 
the hemisphere was measured at each point plotted and 
from all these measurements the average outside radius 
of the hemisphere was obtained. From this value of 
average outside radius and the average thickness of 
the hemisphere obtained as already mentioned, the mean 
radius R g of each specimen was obtained and is shown 
in Table 1. From the drawings of initial and deforma- 
tion profiles mentioned above, radial deflections of 
the deformed specimen were obtained every 10 degrees 
in the meridional direction. These deflection measure- 
ments are shown for each specimen in Table 1-1 to I-l6. 
The detailed measurements taken are not given in this 
report because of their great volume. 

For each specimen a deformation profile in the 
meridional direction was draxm, which was based on the 
radial deflections obtained as mentioned above and on 
the average outside radius of each specimen. 
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The profiles so drawn are shown in figures in the 
end of this appendix and they are drawn from 10 degrees 
up to 170 degrees in the meridional direction. It 
should he mentioned here that due to the large radius 
due to machining in the part of the specimens between 
the hemisphere and the flange (which is shown as R e 
in Table 1) and to the fact that in most specimens the 
thickness of the specimen in this part was found to be 
less than in the rest of the specimen, the deformations 
found there were longer than expected. Prom the profiles 
drawn in the meridional sense the number of peaks was 
measured and is shown in Table 3 . 

In the internally loaded specimens the edge of the 
flange exhibited a downward motion and this had to be 
taken into account in the calculation of the deflections 
of the hemisphere. For this reason the initial measure- 
ments of a point of the flange at a known distance from 
the edge were averaged and this average was subtracted 
from the average of the final measurements of points of 
the flange at the same distance from the edge. In this 
way the vertical movement of the flange at a known dis- 
tance from the edge was found. In order to find the 
vertical motion of the hemisphere due to the deflection 
of the flange, the vertical movement of a point of the 
flange found as described above was multiplied by the 



- 51 - 



ratio of the distance of the cross section of the mid- 
surface of the hemisphere with the flange from the edge 
of the flange to the distance of the measured point 
from the edge of the flange. The vertical motion of 
the hemisphere due to the deflection of the inner edge 
of the flange was subtracted from the vertical deflec- 
tions of all measured points which had a distance from 
the edge of the flange of more than 1.25 in. This dis- 
tance of 1.25 in. was a little less than the calculated 
in each case distance of the cross section of the mid- 
surface of the hemispherical surface with the flange. 

This correction as mentioned above was made to the re- 
sults of all internally loaded specimens except of speci- 
men 7 for which no deflection of the flange was found. 

For the externally loaded specimens measurements 
were taken after deformation in the sense of circles 
initially parallel to the base at two constant dis- 
tances from the edge of the flange. These distances 
(indicated by letter J5 in Figure 1-1) were 1.5 in. and 
2.5 in. for specimen 2 and 1.75 in. and 2.75 in. for all 
other externally loaded specimens. An average value of 
the distance x (as shown in Figure I-l) was obtained 
from the initial measurements at each constant distance 
& and was subtracted from each x distance measured 
after deformation; in this way the horizontal component 
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of the deflection of each measured point was obtained. 
Next the average radius of each circle which was ini- 
tially parallel to the base was obtained by subtracting 
the average x value of initially measured points (at 
constant distance g) as mentioned above from the dis- 
tance of the center of the hemisphere from the edge of 
the flange. This circle was next drawn and the horizon- 
tal deflections of the points measured every 10° were 
plotted on radial lines from the periphery of the circle. 
Thus, 3 6 points were obtained and a deformation profile 
was drawn. For each externally loaded specimen two such 
profiles were drawn corresponding to distances g as 
mentioned above, and they are shown in figures presented 
in the end of this appendix. From these figures the 
number of peaks observed in each profile was counted 
and is shown in Table 3. 

In Figure 1-2 a photograph of specimen 14 after 
deformation is shown. From this photograph the ripples 
due to buckling are clearly shovm, as well as the deep 
crack in the right hand side, which is the reason for 
not measuring the radial deflections of this specimen. 

An explosive thickness of 0.035 in. was used for this 
specimen. 

In Figure 1-3 a photograph of specimen 21 after 
deformation is shown. The ripples due to buckling are 
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also clearly shown here, hut are not as deep as those 
of specimen 14 of Figure II-2. 



- 54 - 



Sample Calculations 



The velocity imparted to the hemispherical speci- 
men by the explosive can be found by dividing the total 
impulse imparted to the specimen by the loaded mass of 
the specimen. The total impulse is obtained by multi- 
plying the specific impulse of the explosive (as obtained 
in Appendix II) times the weight of explosive used. 



where 

1^. = Total impulse 
I g p = Specific impulse 
VJ = Weight of explosive 

Since the angle of the spherical surface 2a which 
is covered by the explosive is known (half of this angle 
is shown in Table 1 as al , the loaded mass of the hemi- 
spherical specimen could be found from the formula: 

p 

Loaded mass = m g = p g 2v (R g ) H(1 - cos a) 

where 

m_ = loaded mass of hemisphere 
s 

p_ = density of hemispherical specimen 

R g = mean radius of hemisphere (given in Table 1) 
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H = average thickness of hemisphere (given in 
Table 1) 

a = half of the angle which is covered by- 
explosive (given in Table 1) 

After the total impulse imparted to the specimen 
and the loaded mass were found as mentioned above, the 
velocity imparted to the specimen V could be found 
as follows: 



v os = V m s 



The nondimens ional impulse parameter X could 
then be found with the existing data. 



where 



X = 



p s^ V os^ R s 



IT 



os 



where 



M = static yield moment of hemispherical 



os 



M 



os 



specimen 

a H 2 
os 

— n — 



o = average static yield stress of the 
o s 

hemispherical specimen. 
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Hence, substituting this value of M QS in equa- 
tion giving X we get: 

. V( V os) 2r s 

X = _ 7~IP 

OS 

For mild steel specimens a QS = 44,967 psi, and 

for aluminum 6C61-T6, a = 40,828; these values are 

os 

averages of the values obtained in tensile tests and 
reported in detail in Appendix III. The densities of 
the materials used are also reported in Appendix III. 

The calculations used to obtain the results for 
test 2 will be shown. For specimen 2 we have the fol- 
lowing data: 

R. = 2.575 inches 
s 

H = 0.07086 inches 
a = 83 degrees 

p s =7.228 x 10" 4 lb f -sec 2 /in. 4 
W = 20.70 grams 

a os = 44 '967 Psi 

Thus, 

m s = 7.228 x 10 -4 lb f -sec 2 /in. 4 x 2 tt( 2.575 in) 2 
x 0.07086 in. x (1 - cos 83 ) 

= 0.0018737 lb f -sec 2 /in. 
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f 



It = I S p xW e = 0.4o69 Tb f -sec/(gram of explosive) 
X 20.7 grains = 8.423 lb f -sec. 

V os = I t //m s = 8 * 42 5 lb f -sec/0. 0018737 lb f -sec 2 /in. 
= 4495.4 in/sec. 

The non-dimensional parameter \ may now be found as 
follows : 

4 x (7.228 x 10' 4 lb f -sec 2 /in 4 ) x (4495.4 in/sec) 2 
\ — i 

44,967 lb f /in 2 x (0.07086 in) 2 
x (2.575 in) 2 = 1.716 

For the calculation of energy ratios E r shown in 
Table 1 the following formula was used: 

2 

o F V 

E - initial kinetic energy _ Ps~~ os 

^r ~ maximum strain' 'energy which " "”2 
can be absorbed elastically a os 

For example, the specimen 2, using a value for young 
modulus of E = 30 x 10^ psi, we get: 

(7.228 x 10 -4 x lb f -sec 2 /in 4 ) x 30 x 10 6 lb f /in 2 
r (44,967 lb f /in 2 ) 2 

x (4495.4 in/sec) 2 = 217 



- 58 - 



For tests 2 6 and 28 for which partial loading of 
the aluminum specimens was used a theory which was deve- 
loped in Reference 9 is used in order to check the ex- 
perimental results of the hemispherical specimen to its 
average thickness as follows: 

W _ v 2 1H 
H ^8 R s a 7 

where 

W = maximum deflection 
H = average thickness of specimen 

p 

v = 1.5(1- cos a) 

a = half of the angle covered by explosive 
X = the same non-dimensional parameter as 
mentioned above 

R g = mean radius of hemisphere 
1 , H 

®7 = T + BIT 

S 



For example, using the data of test 26 we get: 



v = 1.5(1 - cos 2 a) = 1.5(1 - cos 2 45) = 0.75 
X = 767 , R c = 2.5596 in., H = 0.09068 in. 

b 

Q _ 1 1 H_1_l 0.09066 A ca)|)| 

a 7-T + HR~-T + B"x 2.5596 = °' 5044 
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= 0.6311 



W _ 0.75 2 x 767 x 0.09066 in 
TT 48 x 2 . 559 ^ in x 0.5044 



In the same way for test 28 we find: 



-jj- = 1.68 . 



-60- 




- 61 - 




- 62 - 



Figure 



o 
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TABLE 1-1 



SPECIMEN NO. 2 

RADIAL DEFLECTION MEASUREMENTS ON 
SEMICIRCLE 330°-150° 

Angle of Measurement 

From Left Side Measured Radial Deflection 

From Bottom of Flange (Degrees) ( inches) 

10 0.129 

20 0.1275 

30 0.121 

4o 0.129 

50 0.l4o 

60 0.129 

70 0.111 

80 0.118 

90 0.142 

100 0.173 

110 0.176 

120 0.164 

130 0.148 

l4o 0.135 

150 0.124 

160 0.125 

170 0.131 
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Figure 1-5 
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TABLE 1-2 



SPECIMEN No. 4 

RADIAL DEFLECTION MEASUREMENTS ON 
SEMICIRCLE 270° - 90° 

Angle of Measurement 

From Left Side, Measured Radial Deflection 

From Bottom of Flange (Degrees) ( inches 



10 


0.246 


20 


0.193 


30 


0.213 


4o 


0.265 


50 


0.287 


6 o 


0.310 


70 


0.334 


80 


0.335 


90 


0.290 


100 


0.363 


110 


0.373 


120 


0.355 


130 


0.288 


l4o 


0.223 


150 


0.168 


160 


0.130 


170 


0.112 
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Figure 1-7 
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TABLE 1-3 



SPECIMEN NO. 6 

RADIAL DEFLECTION MEASUREMENTS ON 
SEMICIRCLE 270° - 90° 



Angle of Measurement 
From Left Side, Measured 
From Bottom of Flange (Degrees) 


Radial Deflection 
(inches) 


10 


0.178 


20 


0.170 


30 


0.133 


4o 


0.148 


50 


0.179 


60 


0.185 


70 


0.187 


80 


0. 204 


90 


0.222 


100 


0.192 


110 


0.196 


120 


0.178 


130 


0.162 


i4o 


0.135 


150 


0.116 


160 


0.109 


170 


0.120 
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Deformation profiles 

Scale : 5 = f 



Figure 1-9 
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TABLE 1-4 



SPECIMEN NO. 7 



RADIAL DEFLECTION MEASUREMENTS 
SEMICIRCLE 270° - 90° 


ON 


Angle of Measurement 
From Left Side, Measured 
From Bottom of Flange (Degrees) 


Radial Deflection 
( inches) 


10 


o.o4o 


20 


0.028 


30 


0.0251 


4o 


0.026 


50 


0.032 


60 


0.0425 


70 


0.0515 


80 


0.061 


90 


0.071 


100 


0.0625 


110 


0.0535 


120 


0.0475 


130 


0.039 


l4o 


0.034 


150 


0.032 


160 


0.033 


170 


o.o4o 
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TABLE 1-5 



SPECIMEN NO. 9 

RADIAL DEFLECTION MEASUREMENTS ON 
SEMICIRCLE 0° - 180° 

Angle of Measurement 

From Left Side, Measured Radial Deflection 

From Bottom of Flange (Degrees) (inches) 

10 0.117 

20 0.085 

30 0.073 

4 o 0.069 

50 0.080 

60 0.087 

70 0.093 

80 0.097 

90 0 . 108 

100 0.098 

110 0.083 

120 0.071 

130 0.055 

l 4 o 0 . 045 

150 0.0475 

16 o 0.063 

170 0.081 
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TABLE 1-6 



SPECIMEN NO. 11 

RADIAL DEFLECTION MEASUREMENTS ON 



SEMICIRCLE 0 ° - 


180 ° 


Angle of Measurement 
From Left Side, Measured 
From Bottom of Flange (Degrees) 


Radial Deflection 
( inches) 


10 


0.205 


20 


0.1525 


30 


0.1575 


4 o 


0.198 


30 


0.258 


60 


0.276 


70 


0.305 


80 


0.283 


90 


0.273 


100 


0.282 


110 


0.308 


120 


0.287 


130 


0.258 


l 4 o 


0.228 


150 


0.215 


16 0 


0.228 


170 


0.270 
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Figure I -13 
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TABLE 1-7 



SPECIMEN NO. 12 

RADIAL DEFLECTION MEASUREMENTS NO 
SEMICIRCLE 0° - 180° 



Angle of Measurement 
From Left Side, Measured 
From Bottom Flange of (Degrees)_ 


Radial Deflection 
( inches) 


10 


0.060 


20 


0.058 


30 


0.053 


4o 


0.063 


50 


0.077 


60 


0.083 


70 


0.091 


80 


0.087 


90 


0.093 


100 


0.095 


110 


0.090 


120 


0.086 


130 


0.0825 


l4o 


0.068 


150 


0.063 


16 0 


0.081 


170 


0.095 
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Figure I- 1 5 
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TABLE 1-8 



SPECIMEN NO. 13 

RADIAL DEFLECTION MEASUREMENTS ON 
SEMICIRCLE 0° - 180° 

Angle of Measurement 

From Left Side, Measured Radial Deflection 

From Bottom of Flange (Degrees) ( inches) 

10 0.084 

20 0.0725 

30 0.0975 

40 0.132 

50 0.162 

60 0.175 

70 0.1815 

80 0.203 

90 v 0.206 

100 0.189 

110 0.197 

120 0.197 

130 0.184 

l4o 0.1625 

150 0.146 

160 0.150 

170 0.178 
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Figure I— 17 
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Figure I- 18 
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TABLE 1-9 
SPECIMEN NO. 15 

RADIAL DEFLECTION MEASUREMENTS ON 
SEMICIRCLE 0° - 180° 

Angle of Measurement 

From Left Side, Measured Radial Deflection 

From Bottom of Flange (Degrees) ( inches) 

10 0.034 

20 0.0275 

30 0.0278 

40 0.0285 

50 0.0325 

60 0.0328 

70 0.0318 

80 0.0285 

90 0.025 

100 0.0255 

no 0.0152 

120 0.008 

130 0.007 

l4o 0.008 

150 0.0065 

160 0.00 

170 0.012 
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TABLE 1-10 



SPECIMEN NO. 16 

RADIAL DEFLECTION MEASUREMENTS ON 
SEMICIRCLE 270°-90° 

Angle of Measurement 

From Left Side, Measured Radial Deflection 

From Bottom of Flange (Degrees) ( inches) 



10 


0.074 


20 


0.073 


30 


0.066 


4 o 


0.071 


50 


0.074 


60 


0.082 


70 


0.0925 


80 


0.101 


90 


0.108 


100 


0.104 


110 


0.106 


120 


0.093 


150 


0.082 


l 4 o 


0.078 


150 


0.080 


16 0 


0.090 


170 


0.088 



- 89 - 




-9 0 - 



TABLE I -11 



SPECIMEN NO. 17 

RADIAL DEFLECTION MEASUREMENTS ON 
SEMICIRCLE 0° - 180° 



Angle of Measurement 
From Left Side, Measured 
From Bottom of Flange (Degrees) 


Radial Deflection 
( inches) 


10 


0.077 


20 


0.067 


30 


0.059 


4 o 


0.060 


50 


0.069 


60 


0.077 


70 


0.083 


80 


0.080 


90 


0.075 


100 


0.065 


110 


0.058 


120 


0.050 


130 


0.043 


l 4 o 


0.038 


150 


0.037 


160 


0.04-7 


170 


0.0575 
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Figure I - 21 



TABLE 1-12 



SPECIMEN NO. 18 

RADIAL DEFLECTION MEASUREMENTS ON 
SEMICIRCLE 45° - 225° 



Angle of Measurement 
' From Left Side, Measured 
From Bottom of Flange (Degrees) 


Radial Deflection 
(inches) 


10 


0.101 


20 


0.102 


30 


0.092 


40 


0.093 


50 


0.106 


6 0 


0.108 


70 


0.1165 


80 


0.115 


90. 


0.116 


100 


0.112 


110 


0.106 


120 


0.096 


130 


0.088 


140 


0.085 


150 


0.081 


16 0 


0.087 


170 


0.093 
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Figure X-22 



TABLE 1-13 



SPECIMEN NO. 19 

RADIAL DEFLECTION MEASUREMENTS ON 
SEMICIRCLE 0° - 180° 



Angle of Measurement 
From Left Side, Measured 
From Bottom of Flange (Degrees) 


Radial Deflection 
( inches) 


10 


0.0925 


20 


0.058 


30 


0.0565 


4o 


0.071 


50 


0.082 


60 


0.103 


70 


0.112 


80 


0.125 


90 


0.119 


100 


0.097 


110 


0.105 


120 


0.108 


130 


0.104 


l4o 


0.088 


150 


0.078 


160 


0.077 


170 


0.091 
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Figure 1-24 
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TABLE 1-14 



SPECIMEN NO. 21 

RADIAL DEFLECTION MEASUREMENTS ON 
SEMICIRCLE 0° - 180° 



Angle of Measurement 
From Left Side, Measured 
From Bottom of Flange (Degrees) 


Radial Deflection 
( inches) 


10 


0.181 


20 


0.146 


30 


0.156 


4o 


0.202 


50 


0.248 


60 


0.266 


70 


0.291 


80 


0.312 


90 


0.306 


100 


0.298 


110 


0.308 


120 


0.298 


130 


0.276 


l4o 


0.234 


150 


0.204 


-160 


0 . 214 


170 


0.302 
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Figure 1-26 
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TABLE 1-15 



SPECIMEN NO. 26 

RADIAL DEFLECTION MEASUREMENTS ON 
SEMICIRCLE 0° - 180° 

Angle of Measurement 



From Left Side, Measured 
From Bottom of Flange (Degrees) 


Radial Deflection 
(inches) 


10 


0.062 


20 


0.048 


30 


0.042 


4o 


0.067 


50 


0.125 


6o 


0.148 


70 


0.153 


80 


0.1535 


90 


0.1575 


100 


0.142 


110 


0.127 


120 


0.1175 


130 


0.108 


l4o 


0.053 


150 


0.012 


16 0 


0.010 


170 


0.0225 
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TABLE 1-16 



SPECIMEN NO. 28 

RADIAL DEFLECTION MEASUREMENTS ON 
SEMICIRCLE 315 °- 135° 



Angle of Measurements 

From Left Side, Measured ^ Radial Deflection 

From Bottom of Flange (Degrees) ( inches) 



10 


0.072 


20 


0.061 


30 


0.108 


4o 


0.157 


50 


0.156 


60 


0.148 


70 


0.147 


80 


O.I665 


90 


0.175 


100 


0.166 


no 


0.1475 


120 


0.137 


130 


0.133 


i4o 


0.125 


150 


0.088 


160 


0.046 


170 


0.053 
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Figure 1-28 



APPENDIX II 

CALIBRATION OF THE EXPLOSIVE 
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Calibration of the Explosive 



Notation Used 




V 



avg 

x 

Ax 



At 

h 

t 




specific impulse imparted to specimen by the 
explosive 

average velocity of disc for a given test 

displacement of disc from its rest position 
distance between positions of the disc 
time interval between positions of the disc 
nominal thickness of explosive used 
thickness of calibration specimen 
weight of explosive 

weight of disc specimen 

total impulse imparted to specimen by explosive 
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Procedure 



An important part of this work is to determine the 
impulse imparted to the test specimen by the explosive. 
To this end, the velocity of a flat circular plate ac- 
celerated upwards by the explosive was measured. Since 
the mass of the circular plate and that of the explosive 
were known, the specific impulse was calculated as a 
product of the mass of the specimen times the average 
velocity, divided by the mass of the explosive. Hence, 
the units of the specific impulse are dyne-sec per gram 
of explosive. 

The specimens used for the calibration of the ex- 
plosive were steel and aluminum discs at 3 inches dia- 
meter and thicknesses 1/8 inch or 1/4 inch. A Fastax 
(Wollensak WF-2) framing camera (approximately 12,000 
pictures per second) equipped with a standard Fastax 
1,000 cycles per second glow tube was used to record 
the motion of the specimens on 100 feet rolls of East- 
man 4-x Negative Type 7224 film. The set-up used for 
the calibration tests is shown in Figures II-l, II-2 
and II-3 . A 1/4 inch thick foam attenuator (polyeure- 
thene form of polyester type, with a density of 0.032 
gram/cm-^) , cut to the specimen size was attached on the 
specimen with Dupont No. 4684 cement and the explosive 
cut in the same shape was attached on the form using 
the same cement. 
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A steel plate with a 5.5 inches diameter hole in 
its middle was bolted to the table of the blasting cham- 
ber and used as a base. A steel rule with one inch 
intervals marked was erected perpendicular to the plate 
surface at the edge of the hole and was held in place 
by clamps attached on a vertical bar which in turn was 
bolted on the plate. The hole in the plate was just 
large enough for the specimen to go through, thus keep- 
ina a great amount of smoke out of the field of view of 
the camera. In order to hold the disc in position, a 
piece of cardboard with a circular hole in its center of 
diameter slightly smaller than that of the disc was at- 
tached to the steel plate by masking tape. The disc 
rested on this cardboard piece as shown in Figure II-5. 
In all tests made for the calibration the specimens were 
directed upwards. 

The detonation procedure was kept the same as that 
of the tests with the actual specimens. A Detasheet 
leader of dimensions 1/4" x 4o" x 0.015" was folded over 
1/4" at one end and was attached in the center of the 
explosive on the disc. The other end was folded over a 
length of 1.5 inches and then placed in contact with a 
number 6 blasting cap, whose leads were connected to the 
blasting control circuit. 
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The camera was set so that when the circuit was 



actuated it would run for 0.7 seconds in order to reach 
its maximum speed before the blasting cap was detonated. 

A previously calibrated 1,000 cps. timing device recorded 
equally spaced markers on the edge of the film. Two 
1,000 watt spot lights were aimed at the flight path of 
the specimen in order to provide the required illumin- 
ation. All the specimens were painted with white enamel, 
so that their images would be more easily distinguished 
in the film. The speed of the camera in the different 
tests was measured from the films and it was found to 
vary between 11,300 and 14,000 pictures per second. 

The next step was to measure the velocity of the 
specimens using the data taken from the films. Each 
film was first measured to find out the time elapsed be- 
tween two successive frames. This time was based on the 
1,000 cps. timing marks which had previously been checked 
for accuracy with an oscilloscope. The time between 
successive frames was assumed to be constant over the 
length of the film after the explosion. Next each frame 
was projected on a screen placed at such distance from 
the projector so that the images came out in full size. 
Then the images of the specimen as shown in each frame 
were drawn on a sheet of paper which was graduated in 
inches and tenths of an inch and which was attached on 
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the screen. The subdivisions of the measuring rule were 
made to coincide with the subdivisions of the graduated 
sheet., In some tests every image of successive frames 
was drawn while in others every other, every third or 
every fourth image was drawn. A typical graduated sheet 
with uhe specimen images drawn is shown in Figure II-4 
(from test No. 20) . 

Next, the distance of the apparent center of gravity 
of the disc (as shown by its image) from the base was 
measured. The differences in distance from the base be- 
tween each pair of successive images was then found and 
all the differences were averaged. The average velocity 
could then be found by dividing the average distance be- 
tween frames by the time between frames which was already 
calculated from the film. In Figure II-5, a photograph 
of a specimen during a test (test No. 15) is shown; the 
specimen can be seen clearly at the top of the smoke 
which covers the lower part of the graduated steel rule. 

Once the average velocity had been found, the spec- 
ific impulse was found by using the equation: 



Results 



Twenty-six tests were conducted with circular discs 
of which twenty-one gave results and were taken into 
account in the calculation of the average specific im- 
pulse , The average specific impulse is: 



I = 18.172 x 1 (j dyne-sec/(gram of explosive) 
sp 



or 



I *= 0.4o69 

sp 



lh-sec/(gram of explosive) 



Four additional tests were made in order to assess 
the effect of the explosive leader on the specific im- 
pulse. In these tests only the leader was attached on 
the foam attenuator by Dupont No. 4684 cement. In two 
of these tests only an area of 1/4" x 1/4" of the ex- 
plosive way in contact with the foam, while in the other 
two the leader vias folded over a length of 1/4" and then 
the folded part was attached on the foam of the specimen 
(the second case representing the actual way in which 
the leader was attached on the specimen in all the cali- 
bration tests) . In the calculation of the average spec- 
ific impulse mentioned above, the effect of the leader 
was taken into account. 
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As a check with previous calculations of specific 
impulse, in Reference (5), a value of Is = 17.09 x 10^ 
dyne-sec/gram of explosive is reported, in Reference (6), 
a value of 17.69 x 10^ dyne-sec/gram, while in Reference 

2i )± 

(7), values of 19.2 x 10 dyne-sec/grara and 18.42 x 10 
dyne-sec/gram are reported. Thus, it is seen that the 
value of specific impulse obtained herein, lies in the 
range of results obtained previously. 

Table II-l shows the data used in obtaining the 
average velocity in each test. Table II-2 shows the 
data used for the calculation of specific impulse or 
total impulse for all tests made and Table III-3 shows 
the specific impulses of the different tests used in the 
calculation of the average specific impulse. 

In Figure II-6, a graph is presented of the total 
impulse of each test as function of the explosive weight. 
The straight line drawn through the origin has a slope 
equal to the specific impulse referred to above. Thus, 
it is seen that the variation of total impulse with 
explosive weight is linear. In Figure II-7 the cali- 
bration results are shown \tfith specific impulses plotted 
as functions of nominal thickness of the explosive used. 
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Sample Calculations 



For the sample calculations made below the results 
found for test No. 5 were considered. Table II-l shows 
the heights of the successive positions of the specimen 
from the base. It should be mentioned here that these 
heights measurements taken from the motion pictures 
taken during the calibration tests are not exact because 
they neglect the thicknesses of the cardboard (used to 
support the specimens as mentioned above), the attenuator, 
the explosive and the specimen itself. This inaccuracy 
doesn't have any effect on the results however, since 
the average velocity was calculated on the basis of dif- 
ferences in distance between positions of the specimen, 
and not on the basis of distances from the base. The 
average distance between positions of the image of the 
specimen for test No. 5 was found to be 0.377 in as 
shown in Table II-l. The speed of the camera was found 
as 11,500 pictures per second. The positions of the 
specimen were measured every 4 frames and hence the time 

between two positions of the specimen is 

\ 

TT7500 = 57575" seconds 

Therefore, the average velocity in this case will be: 
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V 



0.377 in. 



= 1,084 inches/second 



avg (1/2,875) sec. 



After finding the average velocity, a calculation 
of the specific impulse can he made, since we know also 
the weight of the specimen and the weight of the explo- 
sive which are given in Table II-2. . From Table II-2, 
we have: 



Converting the above found average velocity in centi- 
meter per second we have: 



Finally, with the above data the value of specific 
impulse for test No. 5 is calculated as follows: 



W = 108.35 grams (weight of specimen) 
W g = 1.6 grams (weight of explosive) 



V avg = ^' ( “ ) 84 inches/second = 2,753 cm/sec. 



^p ~ w e 



W D XV 
e 



avg 108.35 grams x 2,753 cm/sec 
l.b grams 




The values for specific impulse calculated in a similar 
way for all tests made are given in the seventh column 
of Table II-2. 



Calculations of Effect of Gravity and Air Drag 



on Specific Impulse 

It is our purpose here to find out how much the 
velocity calculated from the films taken during the cali- 
bration tests must be changed in order to take into ac- 
count the effects of gravity and air drag. 

The following nomenclature is used: 

V c = velocity measured from calibration film 

V = velocity which must be added to V 

g o 

because of gravity effects 
V Q = velocity which must be added to V 
because of air drag effects 
m = mass of calibration specimen 
h = average height of different positions of 
calibration specimen, used in the calcu- 
lation of V Q 

g = acceleration of gravity 
D = air drag 

Considering gravity effects first and equating the total 
energy of the calibration specimen in height zero and h , 
we have: 



“*r ra ( V o + V g ^ 2 = -j- m + mgh . 
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2 

Disregarding the second order term (V ) and solving 

o 

for V , we get: 

o 




Next considering air drag effects alone, and equating 
again the energy at zero height to the energy at height 
h, we get: 

( V o + V a ) 2 = + D * h 

2 

Disregarding again the second order term V and solving 

a 

for V , we get 

ct 




As an example, for test No. 18 with V Q = 4,810 in/sec. 
h = 7.3 in, m = 40.35 grams = 0.0889 lb m we have: 



g - ?a - gft/ aioiaer 12 ln/ft ■ 0-5B6 in/sec 



which corresponds to 0.012^ increase of V Q . For 
the air drag effects we used the following formula: 



D = C D xfv 2 A 
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where 



C D = drag coefficient, which for Reynolds numbers 

-5 

higher than lCr for a disc normal to the 
flow is 1.12 according to reference 8. 
p = air density 

A = area of calibration specimen 

In the case of test No. 18 we have: 

Vd 

Reynolds No. Re = 



where 

V Q = average speed of specimen 
d = diameter of specimen 
v = kinematic viscosity of air 



Re 



It.8ip . in/see x 3 in. _ 5.929 x i 0 5 

I .69 x 10 H ftVsec x 12 in/ft x 12 in/ft 



Hence, C D = 1.12 according to Reference 8. 

p = 0.0683 lt> m /ft 5 

A = = *7p(3)^ = 7.065 in 2 

1.12 x 0.0683 lb m /ft 5 x (-- rg — ft/sec) x 7.065 in 2 
2 x 32.2 lb ft/lb f sec x 144 in 2 /ft 2 

= 9.56 lb f 
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and therefore: 



D>h 9.36 lh f x 7.3 in x 12 in/ft x32.2 lb f' t/lb f sec 
a mV o O.O 889 lb m x 4,810 in/sec 

= 61.743 in/sec 

This corresponds to 1.279^ increase of V Q . In 
test No. 18 and 16, the air drag effects were found 
1.279^ ( as above) and 0.9863 respectively and they 
were taken into account in the calculation of the cor- 
responding specific impulses. In all other cases, 
both gravity and air drag effects were found varying 
from 0.013 to 0.43 and they were neglected. 
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The Effect of Specimen Rotation On Specific Impulse 



The total kinetic energy (KE) of the disc used as 
specimen for the calibration, is equal to the sum of 
its translational energy and the energies of rotation 
about its x, y and z axes, x is the axis lying on 
the plane of the measurement rule, y is the axis 
normal to the plane of the measurement rule and z is 
the axis normal to the plane of the disc and parallel 
to the measurement rule, as shown in Figure II-2. The 
equation for the total kinetic energy is: 




where 

m = mass of calibration specimen 
V Q = average speed of specimen 
I = moment of inertia about x axis 

A 

I = moment of inertia about y axis 
I z = moment of inertia about z axis 
W = angular velocity about x axis 

A 

W = angular velocity about y axis 
«/ 

W = angular velocity about z axis 

In all cases rotation about z axis was neglected 

(W =0). Because of symmetry the moments of inertia 
z 
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about x and y axes are equal (1^ = I = I) . This 
moment of inertia with respect to axes x and y for 
a disc of mass m and radius a is given by: 

2 

-r ma 2 

I = -jj— grams - cm 



if m is given in grams and a in centimeters. 

Hence, the equation for the total kinetic energy becomes 



KE 





(w 2 + w 2 ) 
v x y' 



As an example, let us consider test No. 9. In this 
case, we have: 



m = 107.30 grams, = 5*009 cm/sec h = 7.5 in 
0 =0 (rotation about x axis) 

0 = 22 ° = 0.384 radians (rotation about y axis) 

a = 1.5 in x 2.54 cm/in = 3.81 cm 



I - 2§i = 1 . 07.38^ (3. 81) 2 , 389 . 68 graras _ cm 2 



wf. = 



^ ..2 ^ 0 . 384 x 5.009 ^\ 2 

0 ' w y \ 7V6 ' = 10,195 



sec 



V* 



__ + i^L = 107.38 ( 5 -QO , 9 ) g + 589.68 10 

= (1.5^7 + 0.00198) x 10 9 ergs = 1.34898 x 10 9 ergs 



KE = m 
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If we call the speed calculated after talcing into 

account the rotation effect we will have: 

|-mV 2 = 1.54898 x 10 9 ergs 

V t = ~ = 2 5125 .549 cm 2 /sec 2 

= 5012.5 cm/sec 

Hence we see that the increase in speed due to rotation 
is 0.06984 and is therefore negligible. Similar re- 
sults for rotation effects were also found in the other 
tests and therefore the effects of rotation on specific 
impulse were neglected. 



or 



or 
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Calculation of Effect of Explosive Leader On 



Specific Impulse 



As mentioned above two tests were performed in 
which the leader was folded over a length of 1/4" and 
then attached to the center of the specimen on the at- 
tenuator. These tests are 27 and 28 of Table II-l . 

The total impulse imparted to the specimens was cal- 
culated by multiplying the average speed as shown in 
Table II-l by the weight of the specimen and is shown 

in Table II-2. The average of the two thus obtained 

, 4 

total impulses is 2.405 x 10 dyne-sec. 

In order to find the effect of the leader on the 

specific impulse, the total impulse mentioned above v/as 

subtracted from the total impulses 1^ of all the tests 

and the new total impulse is shown as 1^ in Table II-2 

then these values 1^ so obtained were divided by the 

weight of explosive used for each test and thus the 

corrected value for the specific impulse given as 1^ 

in Table II-2 was obtained. 

As an example, let us consider test No. 1 . The 

initial total impulse before the correction was 
4 

55.72 x 10 dyne-sec as given in Table II-2. The 
corrected total impulse is: 
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1^ = 35.72 x 10 4 - 2.405 x 10 4 = 33.315 x 10 4 dyne-sec 

And since the weight of explosive used in test No. 1 was 
2 grains, the corrected specific impulse is: 

4 . 

i' = — - = 16.657 x 10 dyne -sec/ (gram of 

S p ct 

explosive) 
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TABLE II -1 



DATA USED IN OBTAINING THE AVERAGE VELOCITY OF 
THE SPECIMEN IN EACH TEST MADE 



Test x 

Number (in) 


V 

( Ax) V aVg 

Ax v 'avg At avg (corrected) 

(in) (in) (sec) (in/sec) (in/sec) 



1 4.43 



4.91 


0.48 


5.54 


0.43 


5.77 


0.43 


6.22 


0.45 


6 .65 


0.43 


7.03 


0.38 0.433 3^^ 1,299 1,299 



2 4.29 



4.68 


0.39 


5.10 


0.42 


5.47 


0.37 


5.87 


o.4o 


6.26 


0.39 


6 .68 


0.42 0.3983 2,390 2,390 



3 Film came out very dark, probably due to insuf- 
ficient illumination. 
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TABLE II-l (continued) 



Test 

Number 



4 



5 



x Ax 



3.29 


3.67 


0.38 


4.01 


o.j>b 


4.38 


0.37 


4.73 


0.35 


5.06 


0.33 


5.41 


0.35 


4.08 


4.42 


0.34 


4.80 


0.38 


5.16 


0.36 


5-5^ 


0.38 


5.92 


0.38 


6.30 


0.38 


6.72 


0.42 


3.39 


4.08 


0.69 


4.76 


0.68 


5.41 


0.65 


6.09 


0.68 


6.83 


0.74 


7.54 


0.71 



[ avg 



At 



V 

avR 



V 

avg 

( corrected) 



°- 555 2 > 118 2 > 11Q 



0.377 1,084 1,084 



0.691 575oo 2,075 2,075 
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TABLE II-l (continued) 



Test 

Number x Ax 



4.0 


4.32 


0.32 


4.70 


0.38 


5.07 


0.37 


5.46 


0.39 


5.83 


0.37 


6.20 


0.37 


5.28 


5.57 


0.29 


5.88 


0.31 


6.18 


0.30 


6.46 


0.28 


6.76 


0.30 


7.08 


0.32 


5.38 


6.02 


0.64 


6 .68 


0.66 


7.36 


0.68 


8.04 


0.68 


8.80 


0 .76 


9.50 


0.70 



ve 



0.566 



0.50 



0.686 



At 



V 



avg 



V 

avg 

( corrected) 



1 

57550 



2,068 



2,068 



1 

57W 



1,770 



1,770 
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27B75 



1,972 



1,972 



TABLE II-l (continued) 



Test 

Number x £x 



5.90 


6.44 


0.54 


6.96 


0.52 


7.45 


0.49 


7.98 


0.53 


8.48 


0.50 


8.95 


0.47 


5.78 


6.12 


0.54 


6.50 


0,58 


6.88 


0.58 


7.31 


0.45 


7.77 


0.44 


6.12 


6.54 


0.42 


6.93 


0.39 


7.27 


0.54 


7.57 


0.50 


7.88 


0.51 



(Ax) 



avg 



At 



V - 

V avg 

avg (corrected) 



0.508 57900 2,997 2,997 



o- 59 " srfeo 2 ’* h 2 ' 564 



0.35 2 2,288 2,288 
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TABLE Il-l (continued) 



Test 



Number 


X 


Ax 


13 


4.53 






5.10 


0.57 




5.68 


0.58 




6.20 


0.52 




6.74 


0.54 




7.30 


0.56 


14 


4.89 






5.32 


0.43 




5.76 


0.44 




6.23 


0.47 


15 


6.98 






7.37 


0.39 




7.73 


0.36 




8.12 


0.39 




8.49 


0.38 




8.83 


0.34 




9.20 


0.37 



(“W at 



0.554 



^7750 



0.446 



67750 



0 *^ 16 ?,ooo 



v. 



avg 



3,74o 



3,010 



2,601 



V 

avg 

( corrected) 



3,740 



3,010 



2,601 
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TABLE II-l (continued) 



Test 

Number 



Ax 



(Ax) 



avg 



At 



V 

V , avg X 

avg (corrected) 



16 


5.69 






6.57 


0.88 




7.29 


0.72 




8.08 


0.79 




8.77 


0.69 


17 


Field of 


view of 


18 


4.03 






4.55 


0.50 




4.94 


0.4l 




5.58 


0.44 




5.83 


0.45 




6.30 


0.47 




6.72 


0.42 




7.16 


0.44 


19 


Position 


of spec 


20 


4.58 






5.04 


0.46 




5.52 


0.48 




6.00 


0.48 




6.48 


0.48 




6.99 


0.51 




7.48 


0.49 



0.77 5 -^ 50 - 4,812 4,872 



0.447 



1 

12, ^00 



5,587 5,646 



°- 483 ST15CT 2 '975 
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2,975 



Test 

Number 



TABLE II- 1 (continued) 
(Ax) 



Ax 



avg 



At 



V 

V a v g 

avg (corrected) 



21 


Specimen 


exhib 


22 


6.83 






7.34 


0.51 




7.78 


0.44 




8.25 


0.47 




8.74 


0.49 




9.21 


0.47 


23 


5.97 






6.49 


0.52 




6.99 


0.50 




7.47 


0.48 




7.90 


0.43 




8.32 


0.42 


24 


7.96 






8.90 


1.64 


25 


8.09 






8.67 


0.58 




9.19 


0.52 




9.69 


0.50 



0,476 £7350 3,020 3,020 



°* 47 6735(7 2,985 



1.64 



0.533 



^133 



3,500 



i 

67W 



3,470 



2,985 



3,500 



3,470 



26 Specimen exhibited a large amount of rotation. 
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TABLE II-l (continued) 



Test 



Number 


X 


Ax 


27 


4.45 






4.54 


0.09 




4.63 


0.09 




4.73 


0.10 




4.84 


0.11 




00 


0.09 


28 


5.0 






5.H 


0.11 




5.23 


0.12 




5.33 


0.10 




5.45 


0.12 




5.54 


0.09 


29 


4.4l 






4.48 


0.07 




4.56 


0.08 




4.63 


0.07 




4.70 


0.07 




0.096 



0.108 



0.0725 



At 



V 



avg 



WT 83,2 



wr 93,7 



7 W “ 5 °* 7 



V 

avg 

(corrected) 



83.2 



93.7 



50.7 
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TABLE II-l (continued) 



Test 

Number 



Ax 



(Ax) 



avg 



V avg 

At avg (corrected) 



50 



4.05 

4.18 

4.20 

4.29 

4.58 

4.46 



0.09 

0.08 

0.09 

0.08 

0.08 



0.084 



1 

T5o" 



54.6 



54.6 



Remarks on TABLE II-l 

1. In the above calculations of average velocity, only 
in the tests No. 16 and 18 the correction due to air 
drag was taken into account, since it was found to 
amount to about 1< of the initially calculated values. 
In all other tests the corrections due to air drag 
were negligible. In all tests the corrections due 

to gravity and rotation were found to be negligible. 

2. In Table II-l At refers to the time elasped between 
two successive measurements. 

5. The last four tests are the ones made with the 
leader alone as mentioned above. 
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DATA USED IN COMPUTING THE AVERAGE SPECIFIC IMPULSE 
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TABLE II-2 (continued) 
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TABLE II-2 (continued) 
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TABLE II -3 



CALCULATION OF AVERAGE SPECIFIC IMPULSE 

NO- Isp X IQ* 



1 


16.657 


2 


17.89 


4 / 


19.70 


5 


17.14 


6 


19.18 


7 


18.92 


8 


17.89 


9 


18.55 


10 


18.32 


11 


19.184 


12 


19.06 


13 


17.77 


14 


19.136 


15 


18.10 


16 


17.28 


18 


16.59 


20 


17.10 


22 


19.^2 


23 


18.11 


24 


17.94 


25 


17.68 
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TABLE II-3 (continued) 



The average of the 



Average, specific 



above twenty-one values gives: 
impulse = 18.172x10^ 



Remarks on TABLE II -3 

1. In this table all tests that gave resuls are listed. 

2. The specific impulse 1^ is the corrected specific 
impulse as taken from TABLE II-2. 
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FIGURE II -1 : Set-up for explosive calibration 

inside the blasting chamber 
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FIGURE n -2 : Overall set-up for explosive calibration 



M 
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Masking 

tape 




Figure II- 3 : Bottom view from metal 

baffle of figure IT-1 
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BASE LINE OF MEASUREMENTS 



Figure 11-4: Images of calibration 
specimen from test No. 20 
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Figure 





Figure II- 6 Calibration Results E wc°|ht e 

Wc (grams) 
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APPENDIX -III- 



PROPERTIES OF MATERIALS USED FOR SPECIMEN 

FABRICATION 
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Properties of Materials Used 
For Specimen Fabrication 

Tensile tests were performed in order to find the 
strength characteristics of the materials used in the 
fabrication of the hemispherical test specimens. The 
tensile test specimens were cut from material left after 
the hydroforming of the hemispherical specimens. Since 
the aluminum hemispherical specimens were hydroformed 
in the zero condition of Aluminum 6061, and then heat 
treated to the T-6 condition, the material used in the 
fabrication of the tensile test specimens was also heat 
treated with the hemispherical specimens. The hemi- 
spherical specimens were hydroformed out of metal plates 
of area 1 T x 1» and after the hydroformed hemispheres 
were cut, some of the remaining material was used for 
the fabrication of tensile test specimens. From each of 
the plates used for this purpose two tensile specimens 
were made with fiber direction perpendicular to each 
other; thus the name of the specimens has the same num- 
ber but different letters. Four tensile tests with steel 
specimens and four with aluminum specimens were made. 

The steel specimens come all from the same order of mate- 
rial, while the aluminum specimens come from two dif- 
ferent orders of material. From the tensile tests with 
aluminum specimens only three gave results. 
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An Instron testing machine which belongs to the 
Civil Engineering Department of M.I.T. and which has a 
total load carrying capability of 10,000 lbs. was used 
in conducting all the tensile tests. The machine set- 
ting used for all the steel specimens was 0.1 in/min. 
except for specimen L27 for which a machine setting 
of 0.05 in/min. was used. The machine setting used 
for all the aluminum specimens was 0.05 in/min., except 
for LI 3 for which a machine setting of 0.1 in/in/min. 
was used. The results of all these tests and the speci- 
men characteristics are shown in Table III-l. The yield 
stress for the aluminum specimens was calculated on the 
basis of 0.2^ offset method and is approximate because 
the scale of elongation in the graph produced in the 
Instron machine was not very large. The gauge lengths 
used were about 2 inches for all streel specimens, one 
inch for the aluminum specimens. 

Two typical stress-strain diagrams of the tests 
made are shown in Figure III-l for a steel specimen and 
in Figure III-2 for an aluminum specimen. These were 
obtained from load-elongation diagrams drawn during the 
the tests by the Instron machine, the initial cross sec- 
tional area of the specimen used in order to obtain the 
stress . 
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A chemical analysis for carbon content of the 

steel used for the fabrication of steel hemispherical 

specimens gave a carbon content of 0.21^ . 

Due to time shortage, the densities of the materials 

used were not determined experimentally, but on average 

of the values given in References 3 > 5 and 6 is used. 

—4 p 4 

These average are 7.228 x 10 lb^ sec /in for mild 
steel and 2.49 x 10~^ lb f sec 2 /in^ for aluminum 
606I-T6 . 

Finally, the Young* s modulus used for steel is 

f> f 

30 x 10° psi and that for aluminum 10.5 x 10° psi . 

These values were not measured. 
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SPECIMEN CHARACTERISTICS AND RESULTS OF TENSILE TESTS 
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strain o.30 0.25 0.20 0.15 0.10 0.05 

(in/in) FIGURE EQ - 1 : Tensile test for specimen L19 (mild steel) 



stress (psi) 
5x10 4 




V 
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0.20 0.15 0.10 0.05 

sti ain (in/in) FIGURE jn-2: Tensile test for specimen C13 (aluminum 6061-T6) 



APPENDIX IV 
Explosive Specifications 
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Explosive Specifications 



The explosive used was Du Pont "Detasheet" EL506D. 

Its composition is 63^ pentaerythritol (PETN), 8< nitro- 
cellulose and the remainder an elastomeric hinder. It 
was delivered in rectangular sheets of 10 x 20 inches 
and thicknesses of 0.010" , 0.015" and 0.030" . 

As mentioned in Reference 5> in the present experi- 
mental investigation a lot of problems were experienced 
when explosive thicknesses of 0.010" or 0.020" (which 
is a combination of 2 sheets of 0.010" thickness). In 
cases these thicknesses were used the explosion did not 
propagate all the way across the sheet of explosive used. 
After the tests pieces of the explosive were found intact, 
spread around the blasting chamber, or in other cases the 
explosive on the spherical shell was found burning. How- 
ever, the sheets of explosive with thicknesses 0.015" 
and 0.030" presented no problems at all, as well as two 
sheets of 0.010" and 0.015" put together, the 0.015" 
thickness being on top. 

The Detasheet characteristics according to the manu- 
facturer are given in the following table. 



- 153 - 



Detonation velocity: 


O 

O 

OJ 


meters/sec 


Density 


1.45 


grams /cm-'’ 


Flexibility range 


-30°F 


to 160°F 


Thermal Stability 


25 hours 


250°F 




1 hour 


275°F 




Hot Bar Ignition Temperature 


Instantaneous 


565°f 




30 seconds 


353°F 
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